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Abstract 

We study the classical and quantum cosmology of a universe in which the matter source is a 
massive Dirac spinor field and consider cases where such fields are either free or self-interacting. 
We focus attention on the spatially flat Robertson- Walker cosmology and classify the solutions of 
the Einstcin-Dirac system in the case of zero, negative and positive cosmological constant A. For 
A < 0, these solutions exhibit signature transitions from a Euclidean to a Lorentzian domain. In 
the case of massless spinor fields it is found that signature changing solutions do not exist when 
the field is free while in the case of a self-interacting spinor field such solutions may exist. The 
resulting quantum cosmology and the corresponding Wheeler-DeWitt equation are also studied 
for both free and self interacting spinor fields and closed form expressions for the wavefunction of 
the universe are presented. These solutions suggest a quantization rule for the energy. 
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1 Introdution 



The study of cosmology has always been influenced by the choice of the matter field used to construct 
the energy-momentum tensor in Einstein field equations. The most widely used matter source has 
traditionally been the perfect fluid. However, the ubiquitous scalar field has also been playing an 
increasingly important role in more recent cosmological models as the matter source. This of course 
is not surprising since being a "scalar field" makes it somewhat easy to work with. One may also 
conceivably imagine a universe filled with massless or massive spinor fields as the matter source. Such 
cosmological models have seldom been studied in the literature, and more often than not, in the form 
of general formalisms. In general then, it would be fair to say that cosmologies with spinor fields as 
the matter source are the least studied scenarios. 

A question of interest related to classical and quantum cosmological models is that of signature 
transition which has been attracting attention since the early 1980's. Traditionally, a feature in 
general relativity is that one usually fixes the signature of the space-time metric before trying to 
solve Einstein's field equations. However, there is no a priori reason for doing so and it is now well 
known that the field equations do not demand this property, that is, if one relaxes this condition 
one may find solutions to the field equations which when parameterized suitably, can either have 
Euclidean or Lorentzian signature. 

The notion of signature transition mainly started to appear in the works of Hartle and Hawking 
[1] where they argued that in quantum cosmology amplitudes for gravity should be expressed as the 
sum of all compact Riemannian manifolds whose boundaries are located at the signature changing 
hypersurface. In more recent times, a number of authors have studied this problem when a scalar 
field is coupled to Einstein field equations and shown that the resulting solutions, when properly 
parameterized, exhibit signature transition, see for example [2]. In a similar vain, a classical model 
is studied in [3] in which a self-interacting scalar field is coupled to Einstein's equations with a 
Sinh-Gordon interaction potential. The field equations are then solved exactly for the scale factor 
and scalar field, giving rise to a spatially flat Robertson- Walker cosmology with signature changing 
properties. The case of the non-flat universe is addressed in [4] with a discussion about the conditions 
under which signature transition exists. Also, it is shown in [4] that in the case of a massless scalar 
field this phenomena does not exist. At the quantum cosmology level the same problem is investigated 
in [5] with an analysis pertaining to the exact solutions of the Wheeler-DeWitt equation. Signature 
transition has also been studied in multi-dimensional classical and quantum cosmology in [6] where 
a 4 + d-dimensional space-time is minimally coupled to a scalar field. Finally, it has been used as the 
compactification mechanism in Kaluza-Klein cosmology [7, 8] for a positive and negative cosmological 
constant respectively. 

As a natural extension of the works done in [3] and [4] and as a specific example, we consider a 
spinor field as the matter source interacting with gravity and itself in a Robertson- Walker geometry. 
This is the subject of study in this paper. In general, studying spinor fields coupled to gravity results 
in Einstein-Dirac systems which are not easy to solve. The cosmological solutions of Einstein-Dirac 
systems have been studied in few cases by some authors, see [9] for example and the references therein. 
The fourth reference in [9] is notable in that the quantization of a spinor field coupled to gravity is 
studied in a Robetrson- Walker background. Here we analyze a classical model with spatially flat 
Robertson- Walker cosmology in which a massive spinor field is coupled to gravity. We show that in 
this model signature transition occurs only if the cosmological constant is negative. We also show 
that a massless free spinor field does not lead to signature transition while a massless self-interacting 
spinor field does. 

The quantum cosmology of this model, embodied in the solutions of the corresponding Wheeler- 
DeWitt (WD) equation is studied for both free and self interacting spinor fields. It turns out that the 
WD equation posses exact solutions in terms of Hypergeometric functions. Also, the "zero energy 
condition" imposes a quantization condition on the solutions for the spinor fields. These wavefunctions 
correspond to the classical solutions undergoing signature transition from a Euclidean to a Lorentzian 



domain and, as such, could be useful in understanding the initial condition of the universe. 

The paper is organized as follows: in section 2 we write the Dirac equation in spatially flat 
Robertson- Walker background. Section 3 deals with the energy-momentum tensor of the spinor field 
and the Einstein-Dirac field equations. In section 4 we classify the exact solutions in the cases of 
a free and self-interacting spinor field when the cosmological constant is zero, negative or positive 
respectively and investigate signature transition in each case. The effect of a massless spinor field 
on signature transition is studied in section 5. Section 6 is devoted to the study of the quantum 
cosmology of our model and finally, conclusions are drawn in section 7. 



2 Dirac equation 

Since our aim is to solve the coupled Einstein-Dirac equations in a cosmological background, let us 
concentrate on writing the Dirac equation in a spatially flat Robertson- Walker space-time given by 

ds 2 = -dt 2 + R 2 {t)(dx 2 + dy 2 + dz 2 ). (1) 

As is well known the Dirac equation in curved space-time can be obtained from the Lagrangian 
(c = h = 1) 



^ + r^-^-r (t )^|-7(w), (2) 



where V(ip,ip) is a potential describing the interaction of the spinor field i/j with itself, 7 M are the 
Dirac matrices associated with the spacetime metric satisfying {7^, 7^} = 2g^ u and are the spin 
connections. The Euler-Lagrange equations for ip and then yield 

m + W-^ = 0, (3) 

#S - W + = 0. (4) 

The 7^ matrices are related to the flat Dirac matrices, 7", through the tetrads as follows 

7 M = e^ a , 7m = e^7a- (5) 

For metric (1) the tetrads can be easily obtained from their definitions, that is = e^e^rjaf,, leading 
to 

e« = diag(l, R,R,R), e% = diag(l, 1/R, 1/R, 1/R). (6) 
Also, the spin connections satisfy the relation 

r, = ^ A (^^ + r^) 7 Y. (7) 

Thus, for the line element (1), use of (5) and (6) yeilds 

r = o, r 4 = -| 7 V, (8) 

where 7 and 7* are the Dirac matrices in the Minkowski space and we have adapted the following 
representation 



The preliminary setup is now complete for writing the Dirac equation. From equations (3) and (4) 
we have 
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(10) 



(11) 



Since the metric functions in (1) are functions of t only, we may write ip(r, t) = e ip ' r ifj(t) in (10) where 
P = (Pi>P2,P3) are the separation constants. 

The energy-momentum tensor for the Dirac spinor field as a matter source for gravity can be 
obtained from the definition 

dL 

.'/„,/- (12) 



T — 1- 



which results in 



where 



l in/ 



(13) 



V M V = (fy + r M )V and = - r, 



Since the Einstein tensor associated with metric (1) is diagonal, the energy-momentum tensor of 
its matter source should also be diagonal. With the anzats 4>(r,t) = e tp ' r ip(t) the off-diagonal 
components of become 

T ij = ^{Pi^lj^ + Pj$Jii>) and T 0j = -pj-tpjoif), 

and hence for them to vanish we must have pi = 0, i.e. the spinor field ip depends on t only. Therefore 
equations (10) and (11) simplify to 
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(15) 



The components of the energy-momentum tensor can then be evaluated from equations (13) and (14) 
with the result 
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3 Field equations 

In this section we consider the Einstein-Dirac equation written as 

G^ u = R^ - -Tlg^v + ^Qfiu = KT^ty, ip), (17) 

where Ru, v , 1Z and A are Ricci tensor, scalar curvature and cosmological constant respectively, con- 
structed from torsion-free connections compatible with the metric. The spinor field ip which interacts 
with itself through the potential V(tp,ifj) is a solution of the Dirac equation (14) and is coupled to 
gravity with the energy- momentum tensor given by equation (13). These coupled Einstein-Dirac 
equations must now be solved in a domain that would lead to a spatially flat Robertson- Walker cos- 
mology given by equation (1) with Lorentzian signature (—,+,+,+). However we may parameterize 
the metric in such a way as to allow solutions with continuous transition to a Euclidean domain. To 



this end, we parameterize the metric, as in [3] by adapting the chart {(3, x, y, z] where the hypersur- 
face of signature change would be characterized by (3 = 0. The metric can then be parameterized in 
terms of the scale function R((3) and the lapse function (3 to take the form 



ds 2 = -(3df3 2 + R 2 ((3)(dx 2 + dy 2 + dz 2 ). 



(18) 



It is now clear that the sign of (3 determines the geometry, being Lorentzian if (3 > and Euclidean 
if (3 < 0. For [3 > 0, the traditional cosmic time can be recovered by the substitution t = |/? 3//2 . 
Adapting the chart {i, x, y, z}, we shall write R(t) = R{[3{t)) and also ip{t) = ip{(3{t)). Using equations 
(14) through (18) with units in which k = h = c = 1 one finds the following differential equations 
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(19) 



(20) 



(21) 



(22) 



where a dot represents differentiation with respect to t. We formulate our differential equations in a 
region that does not include (3 = and seek solutions for R and ip that smoothly pass through the 
(3 = hypersurface. The scalar curvature corresponding to our metric is 



(23) 




K = 6 

while the energy density of the spinor field is given by 

p = T 00 = -V$,if>). 



(24) 



It should be noted that equations (19)-(22) are not all independent. We shall see that with a suitable 
choice for the potential V(ip, ip), equation (20) can be obtained by combining equations (19), (21) and 
(22). Also equation (22) is the Dirac equation for the adjoint spinor ip and should not be considered 
as an independent equation. Thus, it is sufficient to solve either equations (19, 21) or equations (20 , 
21) which we shall endeavor to do in the next section after choosing a suitable form for the potential 



4 Solutions 

As mentioned in the last section integrability of the Einstein-Dirac field equations (19)- (22) depends 
on the choice of a suitable form for V(ip,ip). However, this potential should also describe a physical 
self- interacting spinor field. The potential is usually an invariant function constructed from the spinor 
ip and its adjoint Pp. In general, difficulties in treating Einstein-Dirac systems depend directly on the 
form of the spinor ip and the potential V{ip,ip). For simplifying and addressing these difficulties we 
first demand that ip should be a function of t only. This simplifies the Dirac equation (10) into the 
more manageable form (14). We should now choose a suitable self-interacting form for the potential. 
Some of the common forms for V are: V(ip, ip) = mipip representing a free spinor field of mass m, 
V(ip, ip) = mipip J r J^J^ where = ip^^ip and known as the Thirring model, V(ip, ip) = mipip + (ipip) 2 
called Gross-Neveu model and V(ip,ip) = mipip + (iptp) 2 — {iprfip) 2 also known as the chiral Gross- 
Neveu model [10]. Since in this work we are going to study signature transition for the solutions 



of the problem at hand, we concentrate on the two of the simplest forms for V(iJj, ip) which either 
give the exact solutions of the field equations or show signature transition effects. We therefore focus 
attention on the free spinor field and self-interacting spinor field potentials with the latter having a 
term of the form (tpip) 2 . 

4.1 Free spinor field 

This is the case in which the potential has the form V^,^) = rmpip. From equation (16) the 
components of the energy- momentum tensor are Tq = rmf)tp and T\ = 0, and our field equations 
become 

3 (I) - A = -™^, (25) 





(26) 

\ lh I 

^+~V + WV = 0, ( 27 ) 

^ + l^-m^° = 0. (28) 
The solutions of this system can easily be obtained. The results in terms of (3 are listed as follows: 

Case (a) A = 0: in this case integrating equations (25)-(28) results in the solutions 

R(P) = (M/3 3 / 2 + A) 2/3 , (29) 

e ±2/3ira/3 3 / 2 

^ ~ M^ + A (30) 
Also the scalar curvature 1Z and energy density p of the spinor field are found to be 

3M 2 3M 2 , , 

1Z = TT, O = TJ, (31) 

(M/33/ 2 + A) 2 (M/33/2 + A) 2 

where M (being proportional to m) and A are integrating constants. We see that if A ^ 0, the scale 
factor R((3) is not a real function in the Euclidean region (3 < 0. If A = then R{(3) ~ (3 is an 
unbounded function in both Euclidean and Lorentzian regions, passing continuously through (3 = 0. 
However, in this case the functions ip±, 1Z and p all have a singular behavior in (3 = but are well 
defined in the domains (3 > and (3 < 0. The solutions in the case of zero cosmological constant are 
therefore not suitable candidates for exhibiting signature transition behavior. 

Case (b) A < 0: in this case the solutions can be written as 

W) .(il) 1/3 cos2 / 3 (^3 /2 j, (32) 

±2/3im/3 3 / 2 

Y>± ~ -== — — rr , (33) 
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Figure 1: Left, the scale factor for a free spinor field and right, the same figure for a self interacting spinor field, both 
for a typical negative value of A = — 10 -3 with M = g = 1. 
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Figure 2: Left, the curvature scalar for a free spinor field and right, the same figure for a self interacting spinor field, 
both for a typical negative value of A = — 1CT 3 with M = g — 1. 



where M is an integrating constant proportional to m, the mass of the spinor field. The other 
integrating constant is chosen so that R{(3 = 0) = 0. The scale factor R((3) in this case is a real 
function and behaves exponentially for (3 < 0, passing smoothly through (3 = and becomes a 
bounded oscillatory function for (3 > 0. Also, the physical quantities 1Z, ip± and p are regular functions 
both in the Euclidean and Lorentzian domains and pass continuously through the (3 = hypersurface. 
Thus all of the above results show that the solutions in the case of a negative cosmological constant 
exhibit signature transition from a Euclidean to a Lorentzian domain. A plot of the scale factor, 
scalar curvature and energy density in terms of (3 are shown in figure 1 through 3, using equations 
(32) and (34). 

Case (c) A > 0: in this case for integrating the system of equations (25)-(28) we must compare the 

Hubble parameter H = ^ with the value of y^. If H 2 > ^, the solutions of the field equations in 
terms of the parameter (3 read 

m= ^y\^^^ +A y (35 ) 

where M, again proportional to the mass of spinor field m, and A are integrating constants. If we 
take A = 0, the scale factor R{(3) shows regular behavior in (3 < and (3 > regions and also pass 
smoothly through (3 = 0. However, the scalar curvature, energy density and spinor field 
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Figure 3: Left, the energy density for a free spinor field and right, the same figure for a self interacting spinor field, 
both for a typical negative value of A = — 10~ with M = g = 1. 
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all have singularity a at (3 = 0. In the case A / the scale factor (35) is not a well defined real 
function when passing through (3 = 0. If .ff 2 < ^, the above solutions take the form 
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where we have taken the integrating constant A = 0. These solutions resemble the ones for the case 
A < and exhibit signature transition. However, we should note that the energy density of the 
spinor field is negative and thus unphysical. In summery, the above discussion shows that within 
the context of this model, a universe with positive cosmological constant cannot undergo signature 
transition from a Euclidean to a Lorentzian domain through the (3 = hypersurface where its matter 
source is a massive free spinor field. 

4.2 Self-interacting spinor field 

In this section we return to the field equations (19)-(22) with a potential describing a spinor field 
interacting with itself. The form which we have chosen for this potential is that of the Gross-Neveu 
model 

V(ip,ip) = mtpip + X(tpip) 2 , (43) 
where A is a coupling constant. From (16) the components of the energy-momentum tensor read 



T ° = + A(^V) 2 , Tl = Tf = Tf = -\{^f. 



(44) 



The field equations (19)- (22) can now be written as 
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ip + ~ip + (m + 2A^V)7 C V = 0, 



— 3 R - - - 

4> + -— $ - (m + 2AV^)V'7 
2 it 



0. 



If A = 0, the scale factor as a function of (3 obtained from the above equations reads 

R(f3) = (M(3 3 -g) 1 / 3 , 



(45) 

(46) 

(47) 
(48) 

(49) 



where M and g are two constants proportional to the mass of the spinor field and the coupling 
constant A respectively. The other integrating constant is taken to be zero. This function shows 
an unbounded scale factor both in the Euclidean and Lorentzian domains which passes continuously 
through (3 = 0. The scalar curvature and the energy density of the spinor field are found to be 



1 „ , M/3 3 - 3g 
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(50) 



Functionally speaking, these are well defined in both (3 < and (3 > regions, passing smoothly 
through the signature changing hypersurface. However, the question arises as to the validity of these 
solutions since the energy density is again negative in the Euclidean region {(3 < 0), pointing to an 
unphysical matter source. Like the free spinor field case, when A > 0, we are not led to physically well 
defined solutions. In the case of A < the scale factor, being obtained from the system of equations 
(45)-(48) can be written as 
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(51) 



where as before M and g are equal to m and A up to an integrating constant, and the other integrating 
constant is choose so that R((3 = 0) = 0. Also the Ricci scalar and the spinor field energy density are 
obtained as follows 
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A quick look at the functions (51)-(53) shows their regular behavior both in the Euclidean and 
Lorentzian regions without any singularity at (3 = 0. From (51) it is clear that the corresponding 
cosmology is an unbounded universe in the Euclidean domain having a bounded oscillatory behavior 
in the Lorentzian region and passes smoothly through the hypersurface of signature transition at 
(3 = 0. Figure 1 through 3 show plots of the scale factor, scalar curvature and energy density using 
equations (51), (52) and (53). 



5 Massless spinor field 



At this point it would be interesting to take the massless limit of the pervious results. If m = 
and the spinor field is free, i.e. V(ip,ip) = mipip, then from equations (25)-(28) we obtain R(f3) = 
e gC/3 3 / 2 ^ w h ere c i s a constant. This is not a real function for [3 < and thus does not exhibit 
signature transition. In the case of a self-interacting spinor field however, the massless limit can 
be obtained easily by substituting M = in equations (51)-(53), which again results in signature 
changing solutions. To summarize, solutions exhibiting signature transition do not exist when the 
spinor field is free and massless. 



6 Quantum cosmology 

The study of quantum cosmology of the model presented above is the goal we shall pursue in this 
section. For this purpose we construct the Hamiltonian of our model. Let us start with the actin 



<S — J V—9(Lgrav + L matt )d 4 X, 



where 



L 



grav 



n-A, 



(54) 



(55) 



is the Einstein-Hilbert Lagrangian for the gravitational field and L matt represents the Lagrangian of 
the matter source which is given by (2) with g being the determinant of the metric here. Now, by 
substituting (2) and (55) in (54) and integrating over spatial dimensions, we are led to an effective 
Lagrangian in the mini-superspace {R, ip, ip} as follows 



RR 2 + -AR 3 + -R 3 
3 6 



(56) 



Variation of the above Lagrangian with respect to R, ip and ip yields equations (20), (21) and (22) 
respectively. Also, we have the "zero energy condition" given by 



„, dC ■ dC ■ jdC „ 
H = -rR+ -yV + ^ - c = 0, 

dR dip dip 



(57) 



which yields the constraint equation (19). To simplify the Lagrangian (56), consider the change of 
variable u = R 3 / 2 . In terms of this new variable Lagrangian (56) takes the form 



-u + -Ait H — it 
9 3 6 



tpj°ip - ipj°ip - 2V(ip, ip) 



with the corresponding Hamiltonian given by 
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(59) 



Choosing the potential V{ip,ip) = mipip + \(tptp) 2 and noting that equations (47) and (48) allow us to 
write ipip = -fi? = £i and substitute for the terms involving ip and ip, we construct the Wheeler-DeWitt 
equation from Hamiltonian constraint (59) with the result 
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(60) 



where u 2 
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we have taken C = — 1 to obtain a positive energy density for the spinor field, since in the case of a 



16 -X and *(u) is the wave function of the universe. Note that 



positive cosmological constant, our classical solutions (35-42) do not exhibit a physical situation. In 
what follows, we only consider the case where A < 0. 

In order to solve equation (60) we first consider the case g = 0, which is relevant to the free spinor 
field. In this case equation (60) becomes to a Hermite equation and its eigenfunctions can be written 
in terms of Hermite polynomials H n (x) as 

$ n (u) = ^L=e-" u2 / 2 H n (V^u), n = 0,1,2,... (61) 

The "zero energy condition" H = then yields 

27 

m=— (2n+l)w. (62) 

We note that the energy density of the spinor field in this case is p = —mipip = -p-. Thus, in our 
units, m is the total energy of the spinor field. Therefore, equation (62) is the quantization condition 
of the total energy of the spinor field. The general solution of the WD equation in this case can be 
written as 

tf(u) =^c n $ n (u), (63) 

n 

where the prime on the sum indicates summing over all values of n satisfying condition (62). 

For a self-interacting spinor field, that is for g ^ 0, equation (60) after a change of variable v = lou 2 
and transformation tp = v~ l / A (f> becomes 




4 v v' 



+ ^ + -l_^U = 0, (64) 



where k = ^ and p? = + f ■ The above equation is the Whittaker differential equation and its 
solution can be written in terms of confluent Hyper geometric functions M(a, b, x) and U (a, b, x) as 



(j)(v) = e -«/2^+l/2 



dM(n - k + -, 2fi + 1; v) + c 2 U(n - k + -, 2/x + 1; v) 



(65) 



Bearing in the mind that if g = the self-interacting solution should be reduced to solution (61), we 
take c\ = and retain only the functions U(a,b,x) [11]. Therefore, the wave function ty(u) is equal 
to 

*(«) - H) ^ e-^U (i - I + ^±5, 1 + ^-W) • (66) 

A glance at equation (60) shows that the solutions have a behavior of the form ^ ~ e -^« 2 /2 w h en 
u — > oo and ^ ~ ( M 2 )( 1 +v / i+^3)/ 4 i n the limit where u — > 0. Demanding the same limiting behavior 
for solutions (66), one can easily see that the function U(a, b, x) in (66) should reduce to a polynomial 
and this happens when a = —n [11]. This yields a quantization condition for the total energy of the 
self interacting spinor field as follows 

» = f ( 2 „ +1 + ^+l). (67) 

Figure 4 shows ^(-u)! 2 for typical values of the parameters. As it is clear from this figure, the 
wavefunction has a well defined behavior near u = and describes a universe emerging out of nothing 
without any tunneling. Also the largest peak in |^(«)| 2 corresponds to the maximum of the scale 
factor in the classical oscillatory solution as shown in figure 1. 
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Figure 4: Square of the wavefunction for M = g = 1, ui = 0.025 and n = 15. 



7 Conclusions 

In this paper we have studied the properties of the Einstein-Dirac equations in a spatially flat 
Robertson- Walker background, with an eye to the signature changing solutions. We have shown 
that if the cosmological constant is negative, either free or self-interacting spinor fields result in solu- 
tions which admit a degenerate metric in which the scale factor and the physical quantities of interest 
such as scalar curvature and spinor field energy density have a continuous behavior in passing from 
a Euclidean to a Lorentzian domain. This phenomenon also appears when we have a self-interacting 
massless spinor field. The corresponding cosmology in these cases show an unbounded universe in the 
Euclidean region passing smoothly through the signature changing hypersurface and an oscillatory 
behavior in the Lorentzian domain. If the cosmological constant is positive or zero the signature 
changing solutions do not exist. In particular, equations (39-42) show a universe with a positive 
cosmological constant for which although signature transition occurs, the energy density becomes 
negative which does not correspond to ordinary matter. We have also shown that a massless free 
spinor field does not result in such solutions, irrespective of the choice of the potential. 

The quantum cosmology of the model presented above and the ensuing WD equation is amenable 
to exact solutions in terms of Hypergeometric functions. These solutions correspond to the classical 
signature changing solutions and would be of interest in the context of the boundary condition of the 
universe. 
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